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ABSTRACT 

Recent literature has seen an ongoing discussion on the limiting mass of strongly 
magnetized white dwarfs, since such objects may prove to be a source of over-luminous 
type-la supernovae. In an earlier paper, we have presented the mass-radius relation of 
white dwarfs with a strong poloidal magnetic field in Newtonian gravity. The inclu¬ 
sion of effects such as general relativistic gravity and many-body corrections to the 
equation of state can alter the mass-radius relation and the maximum mass. In this 
work we estimate the extent to which these effects may modify the earlier results. We 
find that the general relativistic effects tend to reduce the maximum mass by about 
2% and many-body corrections by another additional ~2%, for an assumed carbon 
composition. We also explore field geometries that are purely toroidal or a mixture of 
poloidal and toroidal and find that the limiting mass of such equilibrium configurations 
can be substantially higher than in the case of a purely poloidal field. 

Key words: stars: white dwarfs — magnetic field — stars: magnetic field — methods: 
numerical — equation of state 


1 INTRODUCTION 

The inclusion of a strong internal magnetic field can sub¬ 
stantially affect the structure and mass-radius relation of 
white dwarfs and raise their limiting mass. Such Super- 
Chandrasekhar mass white dwarfs could be possible pro¬ 
genitors of over-luminous type-la supernovae (Howell et al. 
2006; Hicken et al. 2007; Yamanaka et al. 2009; Scalzo et al. 
2010; Tanaka et al. 2010; Silverman et al. 2011; Tauben- 
berger et al. 2011). Several publications in the recent liter¬ 
ature have addressed different aspects of this issue: for ex¬ 
ample, the effects of magnetic field on the equation of state 
(Das & Mukhopadhyay 2012; Manreza Paret et al. 2015; Zou 
& Meng 2015; Mukhopadhyay et al. 2015), Lorentz force and 
instability (Nityananda & Konar 2014; Coelho et al. 2014), 
possibility of electron capture (Charnel et al. 2013, 2014; 
Vishal & Mukhopadhyay 2014) and total lepton number vio¬ 
lation (Belyaev et al. 2015). In an earlier work (hereafter PI, 
Bera & Bhattacharya (2014)) we have presented the mass- 
radius relation of strongly magnetized white dwarfs with self 
consistent inclusion of Lorentz force. In this work, the grav¬ 
itational field was treated in a Newtonian framework. Sub- 


* E-mail:pbera@iucaa.in 
f E-mail:dipankar@iucaa.in 


sequently, Das & Mukhopadhyay (2015) have pointed out 
that the inclusion of general relativistic (GR) effects may 
substantially reduce the upper mass limit. Here we revisit 
our earlier work with GR treatment of gravity and find that 
the reduction in the upper mass limit is at most ~ 2.2%. 

Another important factor that affects the mass-radius 
relation of the white dwarfs is the equation of state (EoS) of 
the constituent material. In PI we have used a pure Fermi 
degenerate electron EoS to describe the pressure density re¬ 
lation. However, as shown by Salpeter (1961), many body 
effects can have a significant effect on the EoS. In this work, 
we also explore the effects of many body corrections to the 
EoS on the mass-radius relation of magnetic white dwarfs. 
We find that for a white dwarf composed of carbon, this 
leads to a further reduction in the mass limit by about 2%. 

While the above discussion pertains to a purely poloidal 
configuration of the magnetic field, in reality the field con¬ 
figuration can be more complex. We have therefore also ex¬ 
plored configurations with purely toroidal as well as mixed 
toroidal and poloidal fields. Our results in these cases show 
that the mass limit of equilibrium configurations could be 
substantially higher than that for pure poloidal configura¬ 
tion. While such a trend has been reported also by Das & 
Mukhopadhyay (2015), masses obtained by us lie well above 
the limits quoted by them. 
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In our work, we impose a constraint that the maximum 
of the mass density should occur at the center of the star. 
Recent work by Franzon & Schramm (2015), relaxing this 
assumption, finds that even higher mass equilibria may be 
obtained in a doughnut-shaped density distribution within 
the stellar interior. 

In general, a white dwarf may possess fast rotation and 
this may affect its structure. However, we do not include the 
effects of rotation in the present work. The results presented 
here pertain to non-rotating white dwarfs with axisymmetric 
magnetic held distribution. 

The paper is structured as follows. In Section 2 we 
present in brief the method of solving the stellar structure 
equations and describe the equation of state (EoS) used by 
us. The results for different held configurations are presented 
in Section 3. Our conclusions are summarized in Section 4. 


2 EQUILIBRIUM CONFIGURATIONS WITH 
AXISYMMETRIC MAGNETIC FIELD 

To solve for the equilibrium configuration we use the spher¬ 
ical polar coordinates (r, 6 , 0) with origin at the stel¬ 
lar center and the pole located along the axis of symme¬ 
try. We assume that the material has inhnite conductivity, 
satisfying the ideal MHD condition with no electric held. 
The general approach we have followed to solve the stellar 
structure is modeled on the so-called “Self Consistent Field 
method” for both Newtonian and GR cases. These meth¬ 
ods are briehy described in the subsections below. Further 
details on these methods may be found in Hachisu (1986); 
Tomimura & Eriguchi (2005); Lander &, Jones (2009); Bera 
& Bhattacharya (2014); Bonazzola et al. (1993); Bocquet 
et al. (1995); Pili et al. (2014). To obtain the solution in New¬ 
tonian gravity, we use our earlier code (Bera & Bhattacharya 
2014), whereas we use publicly available codes XNS (Buc- 
ciantini X Del Zanna 2011; Pili et al. 2014) and Lorene 
1 to obtain GR configurations. The XNS code hnds self- 
consistently the axisymmetric equilibrium solution of the 
relativistic compact star with rotation and magnetic held 
assuming approximate gravity with a conformally hat condi¬ 
tion (CFC) in GR 3+1 formalism. This code, capable of han¬ 
dling polytropic EoS, has been used to generate equilibrium 
structure of magnetized neutron stars with different held 
geometry (Pili et al. 2014, 2015; Bucciantini et al. 2015). 
It has also been used to study white dwarf configurations 
by approximating the Fermi degenerate EoS as a polytrope 
(Das X Mukhopadhyay 2015; Subramanian X Mukhopad- 
hyay 2015). We have modihed the XNS code to use the full 
form of electron degenerate equation of state, thus avoiding 
the poly tropic approximation. To obtain magnetized white 
dwarf configurations in full general relativistic gravity we 
use the magstar module of the spectral code Lorene which 
iteratively solves for the axisymmetric stellar configuration 
in 3+1 formalism. The output of magstar have been widely 
used to study stellar configurations and as input to the 
study of the evolution of poloidal helds (Liebling et al. 2010; 
Gabler et al. 2012, 2013; Haskell et al. 2014). The advantage 
of using magstar is that it uses GRV2 and GRV3 to test the 


accuracy of the result. GRV2 and GRV3 are similar to the 
virial condition in general relativity in 2d and 3d respectively 
(Bonazzola X Gourgoulhon 1994; Gourgoulhon X Bonazzola 
1994). 


2.1 Newtonian Calculation 


In the presence of a magnetic held, there will be effects of 
Lorentz force on the stellar structure. In the Newtonian limit 
the stellar structure equations are 


-VP = -V$ 9 + - (j X B) (1) 

P P 

V 2 $ 9 = 4t tG P (2) 


where P, p, 4>^, j and B are pressure, mass density, grav¬ 
itational potential, current density and magnetic held re¬ 
spectively. In the ideal MHD condition, the magnetic held 
satishes the following Maxwell’s equations. 

V • B = 0, (3) 

V x B = fjioj. (4) 

where po is the free space permeability. 

In axisymmetry, the set of equations (1 — 4) can be 
transformed, for Fermi degenerate EoS, to the following inte¬ 
gral form using the Hachisu Self-Consistent-Field technique 
(Hachisu 1986; Tomimura X Eriguchi 2005; Lander X Jones 
2009). 


—— E f + $> g = M + C, (5) 

PemB 

where p e is the mean molecular weight per electron, mB is 
the baryon mass, Ef is the Fermi energy, A4 is the function 
which expresses the magnetic influence on the structure and 
C is a constant of integration. The functional form of M can 
be adjusted to give poloidal, toroidal or mixed held config¬ 
urations. Equation (5) is used repeatedly at the grid points, 
along with proper boundary conditions, to converge on the 
stellar configuration iteratively. 

For the poloidal and the mixed held geometry in ax¬ 
isymmetry, Al is a function of the hux-function u — 
r sin OAcf) where A$ is the 0 component of the magnetic vec¬ 
tor potential. For simplicity, we assign to M a functional 
form proportional to u. 


M — const. X u, 
1 


(Vi 


u X 0) + 


( 6 ) 

( 7 ) 


B_ 

r sin 0 

The ^-component of the held is related to an arbitrary 
function of u: B$ — f(u)/rsmQ. Whereas for pure toroidal 
held M and B$ are functions of pr 2 sin 2 0 , and these may 
be written as 


. mK m . 2 • 2 Q\2m — 1 

A4 — — -- 7 (pr sin 0) , 


B^ — K, 


2 m — 1 
(pr 2 sin 2 6) r 


( 8 ) 

( 9 ) 


r sin 0 

here the index m ^ 1 and Km is a constant. 

A stationary solution of the equilibrium structure is ex¬ 
pected to satisfy the stellar virial condition, expressed as 


1 http://www.lorene.obspm.fr/ 


3n + w + je = o 


( 10 ) 
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here 


n: 


contribution of internal energy = 



W : 


gravitational potential energy = 



The projection of the energy-momentum conservation 
equation = 0 on time hyper-surfaces gives a 3-vector 
equation for the matter momentum density. The equation of 
stationary motion for perfect fluid with pure poloidal field 
can be expressed as: 


and 


/ B 2 

— — dV ; 
2/io 


1 dP da _ 1 .(fdAcjy 

e + P dx i dx i e + P^ dx i 


(18) 


equation (18) can be integrated to obtain an expression sim¬ 
ilar to that in equation (5), 


where V is the space volume. 

The stellar configuration obtained from the self consis¬ 
tent field calculation is tested for the virial condition, the 
deviation being expressed in a non-dimensional form, 


_ | 3 n + vk + j(\ 
1vc1 \w\ 


(ii) 


The solution is said to converge when the relative change of 
the configuration parameters fall below a pre-set limit. In all 
the configurations presented by us, the values of \VC\ are 
- 1(T 5 . 


2.2 General Relativistic Calculation 


To obtain the stellar configuration in the general relativistic 
case it is required to solve the Einstein equation 

R^-\g^R=^-T^. ( 12 ) 

Here, the stress-energy tensor T^ v contains matter and elec¬ 
tromagnetic source terms 


= (e + P)vPu v + Pg» v + — 
go 


pRlp" 


-y v F l5 F 


7 <5 


(13) 


where e is the total energy density, u is the fluid 4-velocity 
and F is the electromagnetic field tensor, which can be ex¬ 
pressed using the potential 1-form A as: 


F„ v = - A^„. (14) 

Rfj, u and R are the Ricci tensor and Ricci scalar respec¬ 
tively. The metric tensor components g^ are used to express 
the line element as 

ds 2 = g^dx^ dx u . (15) 


Following 3+1 formalism in axisymmetry the metric tensor 
can be expressed in spherical like coordinates (£, r, 0, (f)) as 

g^dx^dx 1 ' = —N 2 dt 2 


+ A 4 


B 2 r 2 sin 2 0{d(j) — N^dt) 2 + -^-(dr 2 + r 2 d0 2 ) , 

B z 

(16) 


where N(= e a ), N^, A and B are functions of (r, 6) (Boc- 
quet et al. 1995). 

If the non-sphericity of the configuration is small, then 
the metric may be further assumed to be conformally flat 
in spatial dimension and the above expression approximated 
as (Pili et al. 2014), 

g^dx^dx" = -N 2 dt 2 

+ A 4 | dr 2 + r 2 dd 2 + r 2 sin 2 9(dct> — N^dt ) 2 j . 

(17) 


H + a = M + const. (19) 

Here, H is the heat function defined in terms of baryon 
number density n as, 

a( " )= r^TTPwf ( " i)<i " 1 ' (2o » 

where e is the internal energy density. Here too the func¬ 
tional form of M decides the field configuration. 

Equation (19) is used by the XNS and Lorene codes 
to find the stellar configuration. XNS uses the approximate 
line element form expressed in equation (17) whereas the 
full line element form of equation (16) is used by the Lorene 
magstar code. XNS has been designed to give solutions for 
poloidal, toroidal and mixed field conditions but magstar 
handles only a pure poloidal field geometry. Pure toroidal 
configurations using the full metric of equation (16) can be 
similarly obtained (Kiuchi & Yoshida 2008; Frieben & Rez- 
zolla 2012), but we have not attempted this in the present 
work. The functional form of the magnetic term M required 
to model a pure poloidal or a pure toroidal field is similar 
to that mentioned in the Newtonian case. 


2.3 Equation of State (EoS) 


White dwarfs are dense objects and their interior temper¬ 
ature is small compared to the Fermi temperature (Tf) of 
electrons. So the EoS is very similar to that of a zero tem¬ 
perature Fermi electron gas, of which the relation between 
the baryon density (p) and pressure ( P ) may be expressed 
as follows: 


1 3 

x F 


3tt 2 \3 


■ g e m F n e 


h 3 


P = 


3 h 3 


(£f(1 + 2x f )J\ + Xp — sinh 1 x F 
x f {2x f - 3)^T + x^ + 3 sinh 1 x F 


( 21 ) 

( 22 ) 

(23) 

(24) 


Here, n e is the electron number density, A e = x F — 

Pf being the Fermi momentum, e is the internal energy 
density, h is the Plank constant, h— h/2n and c is the speed 
of light. 

Corrections to the above equation of state may arise 
from a number of effects. Coulomb interaction in the ion lat¬ 
tice is one source of such corrections. Salpeter (1961) consid¬ 
ered corrections to the pressure and energy expressions from 
i) Coulomb energy due to ion lattice, ii) Thomas-Fermi ef¬ 
fect due to non-uniform charge distribution, iii) spin-spin 
interaction and exchange energy between electrons and iv) 
correlation energy term. At very high densities the electron 
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Fermi energy becomes high enough to enable neutronisa- 
tion to occur via inverse beta decay. The Neutronisation 
Threshold for He and C are 20.596 and 13.37 MeV respec¬ 
tively (Shapiro V Teukolsky 1983). Beyond this the config¬ 
uration becomes unstable. Hamada V Salpeter (1961) have 
presented the mass radius relation of non-magnetic white 
dwarfs, including these corrections. We investigate in this 
work the effect of these same corrections on the mass-radius 
relation of strongly magnetized white dwarfs. The relativis¬ 
tic Feynman-Metropolis-Teller EoS (Rotondo et al. 2011) is 
an improvement over the Salpeter one, but yields corrections 
very similar to the Salpeter EoS. 

The strong magnetic field may itself modify the equa¬ 
tion of state of white dwarf matter by trapping electrons 
in quantized Landau levels. The effects of this are most se¬ 
vere when electrons inhabit only the lowest Landau levels. 
However, the condition in a realistic stellar configuration 
always distributes electrons over several tens of levels, min¬ 
imizing such impact (Bera & Bhattacharya 2014). Strong 
internal magnetic fields may affect the many-body interac¬ 
tions and modify the Salpeter corrections, however, we have 
not included the effect of Landau quantization and magnetic 
modifications of the many-body interactions in the work pre¬ 
sented here. 


3 RESULTS 

In the following sections, we present the results obtained 
using the methods outlined in the previous section. Ini¬ 
tially different magnetic field configurations from the Newto¬ 
nian calculation are constructed and then we compare pure 
poloidal and toroidal field configurations with their GR re¬ 
sults. Next we present the effects of inclusion of the ion 
interactions in the EoS for poloidal field configurations, in 
general relativistic gravity. We use the baryon mass to rep¬ 
resent the mass of the configuration and adopt a value of 
1.98 85 X 10 30 kg for the solar mass unit (M®). In all GR 
configurations, explored by us, the gravitational mass dif¬ 
fers from its baryon mass by less than 10 -3 . In representing 
the mass-radius relation of non-spherical configurations, we 
use either the polar or the equatorial radius, whichever is 
larger. 

3.1 Magnetic configurations for different field 
geometry : Newtonian calculation 

3.1.1 Poloidal field 

In our calculations, the poloidal field is specified through a 
current distribution which is purely toroidal. The functional 
form of the toroidal current distribution is expressed as 
= rsmOpao. The parameter <a 0 describes the strength of 
the current. An example of a resulting configuration is shown 
in Fig. 1. This object has a central density 2 x 10 13 kg.m -3 
and a central field strength 59.6x 10 9 T. The Newtonian cal¬ 
culation gives a mass of 1.896 M© and a radius of 875.3 km 
for this object. A star with such a strong poloidal field has 
an oblate shape. The ratio of polar to equatorial radius in 
this case is about 0.68 and the ratio of magnetic to gravi¬ 
tational energy is about 13%. Near the equatorial plane in 
deep interior of the star the Lorentz force rises to become 


almost equal and opposite to the gravitational force. How¬ 
ever the plasma beta remains high throughout most of the 
stellar configuration. 

For a given central density and field geometry, a higher 
field strength leads to a larger mass of the equilibrium con¬ 
figuration. The mass-radius relation of strongly magnetized 
white dwarfs is thus dependent on the field parameters (as 
shown in Fig. 4. of PI or Fig. 7). Here we have used the 
central field strength or the ratio of magnetic field energy 
to gravitational energy to represent the field parameter. For 
a fixed central held strength, as the central density is re¬ 
duced the magnetic energy becomes more dominant, caus¬ 
ing the mass-radius relation to deviate strongly from that of 
non-magnetic configurations. We present here solutions in 
which the component of Lorentz force just opposite to grav¬ 
ity never exceeds in magnitude the gravitational force itself. 
Violation of this would move the density maximum away 
from the stellar centre, as shown by Franzon V Schramm 
(2015). 

3.1.2 Toroidal field 

A star with a pure toroidal held (B = is prolate and 

the held is fully conhned within the star with zero surface 
held strength (Fujisawa V Eriguchi 2015). We model such 
a magnetic held distribution by setting m — 1 in equation 
(9), i.e. = KmV sin Op, where iF m is a constant. Pure 
toroidal held configurations have zero held at the center and 
along the axis of symmetry. Therefore in this case we use the 
total magnetic energy as the representative magnetic held 
parameter. A typical stellar configuration arising from such 
a held distribution is shown in Fig. 2. In this example the 
central density is 2xl0 13 kg.m -3 and the ratio of magnetic 
energy to gravitational energy (\M/W\) is about 39%. 

For strong toroidal helds, highly massive stellar config¬ 
urations are obtained as the Lorentz force is able to reduce 
the effective inward gravity for most of the interior (see Fig. 
2 F). The value of plasma beta is less than unity almost 
everywhere except a small region near the axis. Such con¬ 
figurations are likely to be strongly prone to interchange 
instabilities. 

Fig. 3 displays the mass radius relation for white dwarfs 
with pure toroidal held in the interior, as obtained from 
Newtonian SCF calculations. Several curves for different val¬ 
ues of \M/W\ are shown. As the held energy increases, both 
the mass and the radius of the configuration increase in com¬ 
parison to non-magnetic configurations. As seen in the fig¬ 
ure, for \M/W\ = 10% the maximum mass is about 1.76 
M©, while for |A4/W| = 40% the maximum mass reaches a 
value as high as 5.6 M©. 

3.1.3 Mixed field 

In order to investigate structures with a mixture of poloidal 
and toroidal helds, we use a functional form f(u) — — 

Umax ) k+1 , in equation (7). We set k = 0 and M. — aou and 
carry out a Newtonian SCF calculation. Comparing with 
a pure poloidal held structure of the same central density 
and central held strength, we hnd that the introduction of a 
toroidal held increases both the total mass and the oblate¬ 
ness of the conhguration. As listed in Table 3, the mass 
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Figure 1 . The distribution of mass density and magnetic field strength inside a white dwarf configuration with a strong poloidal field 
from the Newtonian SCF calculation. The stellar outer boundary is indicated by a thick line. A. The density distribution shows an 
oblate structure, B. The radial profile of mass density along polar and equatorial directions, C. The 2-d poloidal field distribution, 
D. The radial distribution of magnetic field strength along polar and equatorial directions. E. The Lorentz force component along the 
gravitational force, the dashed white contour marking the region where the Lorentz force is 90% in magnitude of the gravitational force, 
and opposite in direction. F. The distribution of plasma beta, the ratio of matter pressure to magnetic pressure in this configuration, 
the dotted contour indicating the value of 0.5. Plasma beta is seen to be high through most of the interior. 
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Figure 2. Same as Fig. 1 but with pure toroidal field. A. The density distribution shows that the structure is prolate. B. The radial 
density profile along polar and equatorial directions. C. The 2-d toroidal field distribution. D. The equatorial field strength profile. 
E. The Lorentz force component along the gravitational force, the white dashed contour indicating the region where gravity is being 
opposed by a Lorentz force component that is 90% as strong. F. Distribution of plasma beta inside the configuration, the dotted contour 
indicating a value of 0.5. The beta values are much lower in this configuration than in the poloidal case. 
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SCF 

XNS 

Lorene 

CONDITIONS 

grid 

Pc 

Field parameter 

512x512 

2xl0 13 kg.m -3 
Bcore = 59.621 GT 

NR = 500, NTH = 100 
(■ Rmax = 2500) 

2xl0 13 kg.m -3 
k pol = 0.01164, £ = 0 

nr = 129 (inside star) nt = 65 
nr = 65 (outside star) nt = 65 
2xl0 13 kg.m -3 

CFA = 950 

RESULTS 

Mass (M©) 

Reg (KM) 
Rp/Req 

B max (GT) 

\M/W\ 

Virial test 

1.8965 

875.26 

0.6758 

59.621 

0.1288 

|VC| = 7.84x 10 -6 

1.8557 

860.16 

0.6738 

59.893 

0.1305 

1.8540 

865.95 

0.6767 

59.621 

0.1276 

|GRV2| = 2.98X10 -8 

GRV3 = 1.56X10 -8 


Table 1 . Pure poloidal field configurations from Newtonian SCF calculation and from GR codes XNS and Lorene. 



SCF 

XNS 

CONDITIONS 

grid 

1024x1024 

NR = 1000, NTH = 250 
{Rmax = 10000) 

Pc 

2xl0 13 kg.m -3 

2xl0 13 kg.m -3 

Field parameter 

\M/W\ = 0.3928, m = 1 

Km = 7, m = 1 

RESULTS 

Mass (M©) 

5.1611 

4.9490 

Reg (KM) 

7854.6 

7774.6 

Rp/Req 

1.1441 

1.1386 

\M/W\ 

0.3928 

0.3928 

virial condition 

|VC| = 6.27 x 10 -5 

- 


Table 2. Pure toroidal field configurations from SCF Newtonian calculation and from XNS code. 



Figure 3. Mass-radius relation of white dwarfs with strong 
toroidal magnetic field, as derived from Newtonian SCF calcu¬ 
lations. Curves for different values of |A1/VF| are shown. Masses 
here can far exceed the traditional Chandrasekhar limit. 


increases from 1.896 M© to 2.11 M© and then to 2.55 M© 
as the fraction of the field energy in the toroidal component 
is increased from zero to 4.7% and to 7%. The correspond¬ 
ing oblateness also increases, with R p /R eq going from 0.68 
to 0.62 and 0.53 respectively. The field distribution shows 
that the toroidal field is confined near the stellar surface 


with its strength reaching a maximum where the poloidal 
field attains zero magnitude (Fig. 4). The values of plasma 
beta remain high in the interior parts, indicating a relatively 
stable structure. 

3.2 GR effects on white dwarf structure 

In order to quantify the effect of General Relativity on the 
mass and radius of a magnetic white dwarf, we intend to 
compare the configurations obtained using Newtonian and 
GR methods for the same central density and magnetic field 
strength. However this is not straightforward because the 
field strength cannot be directly specified in XNS and Lorene 
codes as input, instead other custom parameters are used. 
The field distribution obtained can only be quantified in the 
output solution. Therefore, an iterative adjustment of input 
parameters to these codes is required to obtain the desired 
interior field strength. We compare here solutions of these 
three codes when their field strengths differ by less than 1%, 
for the same value of central density. 

3.2.1 N on-magnetic configuration: comparison of XNS, 
Lorene and TOV results 

We have modified the GR codes XNS and Lorene to use the 
full form of the Fermi degenerate EoS for the white dwarf 
structure. To validate the modifications of these codes, we 
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Figure 4. Same as Figure 1 but for a mixed poloidal and toroidal field, as in configuration conf. II of Table 3. A. The density distribution 
shows a highly oblate structure. B. The radial distribution of mass density along polar and equatorial directions. C. The field distribution 
in the interior: the poloidal and toroidal components are plotted separately in the left and the right halves respectively. The white contours 
represent the field lines. D. The polar and equatorial distribution of the field strength. E. The Lorentz force opposing gravity rises to a 
maximum of about 75% in this case, the white dashed contour indicating 50%. F. Plasma beta value is very high over most of the star. 
The dotted contour indicating plasma beta value of 0.5 is located close to the stellar surface. 
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Code : SCF 

Mixed field conf. I 

Mixed field conf. II 

CONDITIONS 

Pc 

2xl0 13 kg.m 3 

2xl0 13 kg.m 3 

Field parameter 

Bcore = 59.621, Ko = 5 

Bcore = 59.621 GT, kq = 8 

RESULTS 

Mass (M©) 

2.1128 

2.5552 

Reg (KM) 

913.32 

1000.7 

Rp/Reg 

0.6211 

0.5352 

B max (GT) 

59.621 

59.621 

\M/W\ 

0.1797 

0.2612 

\Mtor/M\ 

4.6753X10- 2 

7.0047X10- 2 

Virial test 

|VC| 2.29x 10 -5 

|VC| = 5.39X10- 5 


Table 3. Mixed field configurations from Newtonian SCF calculation with different toroidal field parameter no for a fixed central mass 
density and central field strength. Here no is a non-dimensional constant that appears in the function f{u). 


TOV 

XNS 

relative difference (XNS & TOV) 

Lorene 

relative difference (Lorene & TOV) 

Mass 1.4158 M© 

1.4155 M© 

0.02% 

1.4167 M© 

0.06% 

Radius 1075 km 

1067 km 

0.74% 

1072 km 

0.3% 


Table 4. Non-magnetic configurations with Fermi degenerate EoS. Here central density = 2xl0 13 kg.m 3 , XNS grid : NR = 500, NTH 
= 100, Lorene grid : nr = 129 nt = 65 (inside star) & nr = 65 nt = 65 (outside star). 



Figure 5. Comparison of Newtonian and GR estimates of masses 
and radii of non-magnetic White Dwarfs. The maximum mass val¬ 
ues are 1.45 Mq and 1.41 Mq in Newtonian and GR calculations 
respectively. The circles represent the percentage of mass modifi¬ 
cation upon inclusion of GR effects. 


first compare the results for non-magnetic spherically sym¬ 
metric configurations with those obtained from the Tolman- 
Oppenheimer-Volkoff (TOV) equation (Tolman 1939; Op- 
penheimer V Volkoff 1939). 

The mass-radius relation of white dwarfs resulting from 
the TOV equation is compared with that from Newtonian 
calculations in Fig. 5. For low mass and larger radius the GR 
mass-radius curve closely follows the Newtonian results, but 
for very compact configurations the curve turns around to 
smaller masses, indicating gravitational instability. To quan¬ 
tify the effect of GR, we compare the masses derived for 
the same central density in Newtonian and GR calculations. 
The relative difference between these values show that the 


change in mass due to GR effects is less than 2% for the 
stable configurations (see Fig. 5). 

The XNS and Lorene results using our modifications 
are compared with the TOV results in Table 4 and Fig 6 A. 
As it is seen, the maximum difference in mass is less than 
0.1% and that in radius is less than 1% for the same cen¬ 
tral density. The density distribution in the stellar interior 
matches closely in all three cases. 


3.2.2 Magnetic configurations in General Relativity 

To illustrate the effect of GR on a magnetic white dwarf 
configuration, we compare in Tables 1 V 2 the results of 
GR and Newtonian calculations for a central density 2x 10 13 
kg.m -3 . The two tables correspond to poloidal and toroidal 
field geometries respectively. The poloidal field configuration 
is constructed by specifying a toroidal current distribution 
of the form — (e + P)k po i. We set k po i = 0.01164 and 
£ = 0 in XNS to yield a central field strength of 59.6 GT. 
This combination of central density and field strength yields 
a near-limiting mass configuration. In Lorene code the same 
field strength is achieved by setting the parameter CFA = 
950. With these settings the GR calculations yield masses 
1.855 M© (XNS) and 1.854 M© (Lorene), with radii of 860 
km and 866 km respectively (Table 1). Comparison with 
the corresponding Newtonian configuration shows that the 
mass of this configuration is reduced by ~2% when general 
relativistic effects are included. The density and field distri¬ 
butions of the GR and Newtonian configurations differ very 
little, as seen in Fig. 6B. 

We construct a number of poloidal field configurations 
with different central densities and field strengths to evalu¬ 
ate the difference in mass derived from Newtonian and GR 
treatments. The results are shown in Fig. 7. As expected, the 
difference is larger for more compact configurations. The dif- 
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Figure 6. A.: Comparison of the density profiles obtained from 
TOV, Lorene and XNS calculations of the non-magnetic white 
dwarf configuration presented in Table 4. The difference in the 
density distribution with respect to the TOV solution is displayed 
here. The XNS calculation was terminated at a lower density 
threshold of 10 7 kg.m -3 , about 10 -6 times the central density. 
The peak discrepancy between the density profiles is less than 
1%, while the total mass differs by less than 0.1%, validating 
the modifications we have made to the Lorene and XNS codes. 
B.: The difference between the density (upper panel) and field 
strength (lower panel) distributions obtained from GR and New¬ 
tonian calculations of a white dwarf configuration with strong 
poloidal magnetic field (Table 1). XNS and Lorene results, and 
polar (p) and equatorial (eq) profiles are shown separately as in¬ 
dicated. Deviations in the profiles due to GR are limited to within 
~ 1% of the maximum, while the total mass differs by ~ 2%. C.: 
Same as panel B but for pure toroidal field (Table 2). The GR 
correction here goes up to 10% of the maximum in the profiles 
and ~ 4% in total mass. 
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Figure 7. Mass-radius relation of white dwarfs with poloidal 
magnetic field. Lines are derived from Newtonian SCF calcula¬ 
tions and symbols along them indicate the assumed central field 
strength (expressed in units of B c = 4.41 xlO 9 T). Configura¬ 
tions with a fixed central magnetic field start to deviate from 
the non-magnetic curve as the central density is reduced. Con¬ 
figurations with a fixed ratio of magnetic energy to gravitational 
energy are described by a mass-radius relation similar to that of 
non-magnetic ones, but shifted to a higher mass. The curve for 
|A4/Wj = 10% is shown, it has a maximum mass of about 1.79 
Mq. The squares indicate the percentage of modification in mass 
if general relativistic effects are included. 


ference is large also for stronger magnetic fields as they yield 
more massive configurations. 

For a pure toroidal field with |.A4/VF| ~ 39% and the 
central density the same as above, the resulting XNS config¬ 
uration has a mass of 4.95 M©, which is ~ 4% smaller than 
the Newtonian result (Table 2). The radial distribution of 
density and field agree well in the two cases (Fig. 6C). 

We have not carried out a similar comparison of GR and 
Newtonian results for mixed field confgurations, but expect 
the general trend of reduction in mass by a few percent to 
apply in this case as well. 


3.3 Salpeter EoS 

So far we have used ideal Fermi degenerate EoS to obtain 
configurations with different field geometry. We now investi¬ 
gate the consequence of including many body effects in the 
EoS, a la Salpeter (1961). Table 5 presents a comparison of 
the pure poloidal field configurations using ideal Fermi EoS, 
and those using Salpeter EoS for Helium and Carbon compo¬ 
sitions respectively. The computations here are carried out 
in full GR, using the Lorene code. All three configurations 
have the same central density, 2xl0 13 kg.m -3 and the cen¬ 
tral magnetic field strength is adjusted to be as close to each 
other as possible (~ 59.6 GT). We see that the mass of the 
configuration is reduced upon inclusion of the many-body 
effects. The reduction is ~ 0.7% in case of He composition 
and ~ 2% in case of Carbon. In the absence of the mag¬ 
netic field, the same central density would give masses of 
1.403 M© and 1.386 M© with Salpeter EoS, for He and C 
composition respectively. The central density chosen is be¬ 
low the neutronisation threshold in either case. However if 
the composition is Fe, then the neutronisation threshold is 
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1.1 x 10 12 kg.m 3 which yields a maximum stellar mass of 
1.45 M 0 with strong poloidal field. 


4 CONCLUSIONS 

In this paper, we have computed strongly magnetized ax- 
isymmetric white dwarf configurations, considering poloidal, 
toroidal and mixed poloidal-toroidal field geometries. In all 
cases, the increase of magnetic energy in a configuration in¬ 
creases its mass. We carry out our calculations in both New¬ 
tonian and GR, and quantitatively evaluate the magnitude 
of corrections introduced by GR effects. For a typical white 
dwarf, the radius is about a thousand times larger than its 
Schwarzschild radius, indicating that the GR corrections are 
expected to be small. From our calculations we find that the 
correction to the mass of the configuration is at most 2-4% 
depending on the field geometry. Our main results are : 

(i) The Lorentz force has a strong effect on the structure 
of the stellar configuration. The magnitude and the orienta¬ 
tion of the Lorentz force depend on the field geometry. For 
pure toroidal field configuration the magnetic effects could 
be extremely large and configurations with mass more than 

5 M 0 can be obtained. 

(ii) The inclusion of general relativistic effects in the cal¬ 
culation of stellar structure, in general, reduces the stellar 
mass and radius. The percentage of modification is depen¬ 
dent on the compactness of the configuration. For a pure 
poloidal field geometry, GR reduces the maximum mass by 
about 2% and for a pure toroidal field geometry ithe maxi¬ 
mum correction is about 4%. 

(iii) The inclusion of many body interaction effects in the 
EoS also modifies the stellar configuration and reduces the 
maximum mass. The extent of reduction in mass is depen¬ 
dent on the composition of the white dwarf. For Carbon or 
heavier atoms, this correction is larger than that due to GR 
effects. 

The configurations presented here satisfy conditions of 
equilibrium, but are not tested for stability. We restrict 
the interior density below the neutronization threshold and 
the GR instability threshold. However the strong magnetic 
field in the interior may itself drive interchamge instabilities 
in short time scales (Bera & Bhattacharya 2014). Further, 
pure poloidal and pure toroidal magnetic field distributions 
are not stable by themselves and may evolve over Alfven 
time scale (Markey & Tayler 1973; Tayler 1973; Mitchell 
et al. 2015). An appropriate mixture of poloidal and toroidal 
fields may provide stability (Braithwaite 2009; Braithwaite 

6 Nordlund 2006). The mixed field configurations explored 
by us have, however, a relatively small toroidal compo¬ 
nent (\Mtor/M\ ^ 7%), which may not be sufficient to 
achieve stability. It appears unlikely that mixed field con¬ 
figurations with 10%) magnetic energy in the toroidal 
component could be constructed in axisymmetry (Cioffi & 
Rezzolla 2013; Armaza et al. 2015), and more complex, non- 
axisymmetric structures may be responsible for the ^ 2 M 0 
white dwarf progenitors of over-luminous Type la super¬ 
novae. 
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Code : Lorene 

Ideal degenerate EoS 

Salpeter EoS for He 

Salpeter EoS for C 

CONDITIONS 

Pc 

2xl0 13 kg.m 3 

2xl0 13 kg.m 3 

2xl0 13 kg.m 3 

Field parameter 

CFA = 950 

CFA 954 

CFA = 959.5 

RESULTS 

Mass (Mq) 

1.8539 

1.8405 

1.8194 

Reg (KM) 

865.95 

861.47 

856.28 

Rp/Reg 

0.6767 

0.6752 

0.6749 

B max (GT) 

59.621 

59.607 

59.592 

\M/W\ 

0.1276 

0.1284 

0.1285 

Virial test 

GRV2 = 2.98xl0 -8 
GRV3 = 1.56X10- 8 

|GRV2| = 2.694xl0 -7 
GRV3 = 1.696xl0 -7 

|GRV2| = 2.336x 10 -7 
GRV3 = 2.135X10- 7 


Table 5. Pure poloidal field Configurations using ideal degenerate EoS and Salpeter EoS for Helium and Carbon from relativistic code 
Lorene. 
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